Directed Fixed Energy Sandpile Model 
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We numerically study the directed version of the fixed energy sandpile. On a closed square lattice, 
the dynamical evolution of a fixed density of sand grains is studied. The activity of the system shows 
a continuous phase transition around a critical density. While the deterministic version has the set 
of nontrivial exponents, the stochastic model is characterized by mean field like exponents. 
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Spontaneous emergence of long ranged spatio-temporal 
correlations under a self-organizing dynamics, in absence 
of a fine tuning parameter is the basic idea of Self- 
organized Criticality (SOC). Sandpile models are the pro- 
totype models ofSOCllllHIlHlSSIEI!- In S&nd " 
pile models an integer height variable hi representing the 

number of grains in the sand column is associated with 
every site of a regular lattice. The system is driven by 
adding unit grains of sand at a time. When the height 
hi > h c the sand column topples and it looses some 
grains which are distributed among the neighboring sites 
This creates an avalanche of sand column topplings 
and the extent of such activity measures the size of the 
avalanche. Sand grains go out of the system through a 
boundary so that in the steady state the fluxes of inflow 
and outflow currents balance. 

On the other hand a fixed energy sandpile (FES) 
(ill ITU H^ . Hi^ is a sandpile model within a closed sys- 
tem. Therefore the total mass of sand in this system is a 
conserved quantity. The control parameter is the density 
C of grains. A stable state has heights varying from to 
h c — l and is called an inactive state, whereas, any height 
configuration that has at least one unstable site is said 
to be in the active state. The dynamics of the system 
starts with a random distribution of N = C,L 2 grains. 
Initially some sites may be unstable, which topple. Con- 
sequently some of the neighboring sites may topple again 
and the activity continues. For an infinitely large system 
there exists a critical threshold Cc such that if C < Cc the 
activity terminates and the system gets absorbed in an 
inactive state where as for C > Cc the activity of the sys- 
tem fluctuates but maintains a steady mean value |10| . 
Therefore Cc is the critical point of a continuous phase 
transition from an absorbed phase to an active phase. 

After the dynamics starts, the system takes some time 
to relax to the steady state. The activity at a certain 
time is measured by the fraction p of lattice sites which 
are unstable at that time. In general the mean activity 
(p) is a function of the density C i-e., the deviation from 
the critical point A = C — Cc and also the system size L. 
The simultaneous dependence of activity on A and L is 
expressed by the following scaling form: 
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where Q{x) is an universal scaling function such that 
Q{x) — * x 13 when x » 1. This implies that for a cer- 
tain range of A if L is so large that L 1 / 1/± A >> 1 then 
(p(A, L)) is independent of L and depends solely on A 
as (p(A,L)) ~ A' 3 . Naturally (3 is the order parameter 
exponent for the transition. On the other hand when 
x « 1, Q{x) — > constant, independent of both A and L 
implies that right at the critical point C = Cc the order pa- 
rameter varies with the system size as: (p(L)) ~ L~P/ U± , 
independent of A. 

To which universality class the FES model exponents 
should correspond to? Intensive research has been done 
to study the universality class of the phase transitions in 
FES models. It has been suggested that FES belongs to 
the universality class of the linear interface models (LIM) 
but not of that of the Directed percolation (DP) univer- 
sality class [lfj, [ill, [13, LU1 • DP is generic for continuous 
absorbing state transitions in the absence of a conserva- 
tion law where as in FES there exists a conserved field 
which is the density and it couples the order parameter, 
i.e., the mean activity. We, in the present study like to 
examine if an explicit application of the directional bias 
to the FES system makes the system behave as DP or it 
results to another new universality class. 

Application of a global directional bias onto a system 
has been proved to have strong effect on the critical be- 
haviors of various models in Statistical Physics. In di- 
rected systems, degrees of freedom of the individual ele- 
ments is reduced, which shrinks the configuration space 
of the system compared to undirected system. As a re- 
sult a directed system is simpler and easily tractable an- 
alytically. Examples include Directed percolation [T^ |. 
Directed Sandpile Model [HE[l3, Direct ed River net- 
works 0| and Directed Self-avoiding walks etc. 

We study here the directed fixed energy sandpile 
(DFES) models on an oriented square lattice placed on 
the x — y plane, periodic boundary conditions are im- 
posed along both the directions. A preferred direction is 
imposed onto this system along the — y direction. The 
critical height of stability of a sand column is fixed at 
h c = 1, and on a toppling two grains of sand are dis- 
tributed to the two neighboring sites along the preferred 
direction i.e., at the lower-left (LL) and lower-right (LR) 
positions. 
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FIG. 1: Filled circles denote the sites on the toppling front in 
an infinite avalanche of DDFES on an oriented square lattice 
of size L — 32. This avalanche is periodic and has a period 
64. Horizontal empty circles denote positions of the TF in 
another 63 time units. The two end sites of the TF fluctuates 
but maintains a mean distance of L/2. 

Two different versions of the model are studied accord- 
ing to the rule by which the two grains are distributed in 
a toppling: (i) both the LL and LR sites get one grain 
each, which we call as the deterministic directed fixed en- 
ergy sandpile model (DDFES), (ii) each of the two grains 
is distributed randomly to any of the LL or LR positions, 
and this version is called the stochastic directed fixed en- 
ergy sandpile (SDFES) model. 

We shall first discuss the properties of DDFES. The 
critical point can be arrived at from an inactive state 
by adding grains one by one on an initial empty lattice 
followed by the relaxation of the avalanche. On the av- 
erage both the size and the life times of the avalanches 
increase as the density grows. The toppling front (TF) of 
an avalanche is a set of horizontal contiguous sites which 
travels downward with unit speed. The length of the TF 
however fluctuates (Fig. 1). If at an intermediate time 
the TF has n sites, then at the next time step its length 
can be only n— 1, n or n+1. For a finite avalanche the TF 
first grows from a single site to a certain length and then 
shrinks to zero. The set of sites covered by the left and 
right end sites of TF are the paths of two annihilating 
random walkers ^5|. For a finite avalanche they meet 
and annihilate, however for an infinite avalanche these 
two random walkers cannot meet and the best possible 
way it can be ensured if they can maintain a distance of 
L/2 on the average from each other. 

Therefore the minimum possible sustained activity for 
DDFES in a system of size L is p = 1/(2L). What is 
interesting is, on such a system if the density is slightly 
increased the avalanches created by the additional grains 
die away and the system maintains the activity of the 
infinite avalanche. In a sense the system gets locked with 
this activity for a certain range of grain density. However 



FIG. 2: The mean activity (p) in a system of size L — 64 
which grows with the density in a step like manner where the 
step heights are Ap = 1/(2L). 

on increasing the density even further, a second infinite 
avalanche is created and both the TFs run simultaneously 
resulting a sudden jump in the activity by doubling its 
magnitude to 1 / L. This continues for some range of grain 
density which ends at another jump in activity to 3/ (2L). 
Thus in general the variation of activity is discrete and 
has a step like variation with step heights 1/(2L) (Fig. 
2). As the system size increases the step height decreases 
to zero and the variation of p with £ becomes more and 
more smooth. Similar step-like behavior of the order pa- 
rameter was also observed in [2fJ. 

The configurations at the steady state are periodic and 
the same detailed distribution of grain numbers at all 
sites repeat at regular intervals of time. The periodic 
time is always multiples of L in a L X L system. For 
small systems this period has different values for different 
initial configurations but in most (about 95%) cases the 
period is 2L and rarely L, 3L, AL etc. However, for bigger 
system sizes e.g., for L = 512, 1024 or more the period is 
always 2L. This helps to calculate the order parameter. 
Given an initial distribution of grains, it therefore needs 
to find out the mean activity over only a period and then 
average over many initial configurations. 

The critical density £ c actually has a system size de- 
pendence. To study this variation we start with a closed 
empty lattice of size Lx L and go on adding sand grains 
one by one at randomly selected lattice sites similar to 
what is done in an open sandpile. The dynamics of the 
avalanche is followed for each sand grain added. The 
mean avalanche size increases with the density of grains 
in the system and at a certain £ = Cc{L) depending on the 
sequence of randomly selected sites at which the grains 
were dropped, the activity does not stop any more and 
an "infinite" avalanche continues for ever. In practice, in 
our simulation we followed an avalanche up to a certain 
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FIG. 3: The variation of the deviation of the critical density 
Cc(L) of a system of size L from its value £ c at the infinitely 
large system is plotted with L _1/V H . For DDFES we obtain 
C c » 0.4115,^1, = 1/0.361 « 2.77. 

relaxation time T = 10 6 for L < 128 and T = 5 x 10 6 
for L > 128 to declare the avalanche as infinite. Repeat- 
ing this simulation a large number of times, every time 
starting from an empty system, we calculate the average 
critical density (£ C (L)). These values are then extrap- 
olated as: (Cc(L)) = Cc + AL~ X I V W as shown in Fig. 3 
to obtain ( c = 0.4115 ± 0.002 and l/v\\ = 0.361 giving 
V\\ w 2.77. We also esimated Q c by the scaling plot of 
mean avalanche size (s(L))L~ - 05 vs. AL 003 which is 
also consistent with our estimate of ( c — 0.4115. 

A system of size L is filled initially with density £ by 
randomly distributing QL 2 grains. We allow the system 
to evolve up to a relaxation time T after which the ac- 
tivity is measured at every time step. To measure the 
mean activity (p(A, L)) for a slightly higher density (+S( 
we take the advantage of the fact that the system dy- 
namics is deterministic. On the same initial distribution 
of grains corresponding to the density £ another (5C,)L 2 
grains are randomly added. This ensures that if certain 
density gives sustained activity, its higher density nec- 
essarily gives a non-stop activity. These measurements 
are then repeated for different system sizes. In Fig. 4 
we show the scaling of the order parameter on double 
logarithmic scale. Plotting (p(A, L))L with L°- 55 A we 
observe a nice data collapse for system sizes L =128, 256 
and 512. Comparing with the Eqn. 1 we conclude /3/v± 
= 1 and 1/i/_l=0.55. This implies that j3 = u x « 1.82. 

The analysis so far enables us to estimate the dy- 
namical exponent z — V\\/v_\_ « 1.52. This value of 
the dynamical exponent is directly verified by measur- 
ing the survival probability. The survival probability 
P(t) that the initial activity in a random distribution 
of grains survives a time t has an exponential distribu- 
tion as: P(t) ~ exp(—t/r). At the critical point ( c the 
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FIG. 4: Scaling of the order parameter (p(A,L)} with the 
deviation A = ( — £ c from the critical point. From this data 
collapse and Eqn. (1) we find, f3/v± = 1 and l/v± — 0.55. 



characteristic time is a function of only the system size as: 
t{L) ~ L z where, z is the dynamical exponent z = V\\/v\_ 
of the system. Therefore we calculate the average sur- 
vival time (t c (L)) which is also proportional to L z at Q 
for different system sizes. This is done again by dropping 
grains of sand one by one into a closed system and cal- 
culating the life time of the largest avalanche before the 
system gets locked into an infinite avalanche. Averag- 
ing over many initial configurations the largest life-time 
(t c (L)) is plotted in Fig. 5 for on a double logarithmic 
scale. The slope of the straight line gives the value for 
the dynamical exponent z — 1.49 ±0.05 for DDFES com- 
pared to 1.52 obtained previously. 

In the stochastic directed fixed energy sandpile model 
the critical density £ c is found to be very close to 0.5. 
The order parameter has a highly linear variation with 
A as: p(A) = AA where A w 0.46. A similar to DDFES 
calculation of the system size dependent critical density 
after extrapolation (£ C (L)) = Cc + A'L -1 /"!! gives uu w 1 
and A' « 2.4. The plot of (p(A,L)) vs. A is a very nice 
straight line and the plot of data for different system sizes 
fall on top of one another. We conclude that j3 w 1 and 
v± ~ oo. 

The roughneing of the associated interface in our FES 
models is studied. If Hi(L, t) denotes the number of top- 
plings upto the time t then the set of Hi(L,t)s for all 
i represents an interface. For DDFES the width of this 
interface fluctuates periodically but its average grows as: 
W(L,t) - L a Q(t/L z ) where we find for a = 0.31 and 
z = 1.6 in comparison to the Linear Interface Model re- 
sults a = 0.75 and z = 1.56 t 22|. 

Finally we study the DFES models on oriented square 
lattices of rectangular shapes, the longer sides being par- 
allel to the preferred direction. For DDFES model the 
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FIG. 5: The average of maximal life-time (t c (L)) of the 
avalanche prior to the infinite avalanche is plotted with L. 
The slope gives a measure of the dynamical exponent z = 
1.49. 



FIG. 6: Two snapshots of height configurations in the di- 
rected fixed energy sandpile model on a 32 x 64 oriented 
square lattice, downward direction being the preferred direc- 
tion. Stochastic DFES is shown on the left where as the 
deterministic DFES is shown on the right. Active sites are 
shown by filled circles, open circles denote sites with height 1 
and vacant sites are not indicated. 



Model 




P 


v±_ 


"II 


DDFES 


0.4115 


1.82 


1.82 


2.77 


SDFES 


0.5 


1.00 


oo 


1.00 


BTW FES 


2.125 


0.7 


0.90 


1.49 


Manna FES 


0.71695 


0.64 


0.82 


1.29 


DP 




0.583 


0.733 


1.295 



TABLE I: Comparison of critical points and exponents for 
different models of fixed energy sandpiles. Exponent values 
for BTW and Manna sandpiles are taken from |12|. DP expo- 
nents from |2l|. 

TFs are contiguous sites covering the transverse direction 
in the form of rings. These toppling rings are perfectly 
stable, once formed they never change in shape. As den- 
sity increases, the number of such rings increases. On the 
other hand for SDFES, the toppling sites are randomly 
scattered throughout the system (Fig. 6). 

To summarize, we have studied the directed version of 
the fixed energy sandpile on the oriented square lattices. 
Like isotropic FES, our directed FES also shows a con- 
tinuous phase transition from an absorbed phase to an 
active phase. Two versions of the model are studied. In 
the deterministic FES, the grain number configurations 
are periodic and repeats at regular time interval of 2L. 
For this model the critical points as well as the criti- 
cal exponents are found to be non-trivial and belong to a 
new universality class. The other version has the stochas- 
tic toppling dynamical rules and exponents of mean-field 
nature are found for this model. 
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